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Abstract 

We have studied the ordering kinetics of a two-dimensional anisotropic Swift-Hohenberg (SH) model 
numerically. The defect structure for this model is simpler than for the isotropic SH model. One finds 
only dislocations in the aligned ordering striped system. The motion of these point defects is strongly 
influenced by the anisotropic nature of the system. We developed accurate numerical methods for following 
the trajectories of dislocations. This allows us to carry out a detailed statistical analysis of the dynamics of 
the dislocations. The average speeds for the motion of the dislocations in the two orthogonal directions obey 
power laws in time with different amplitudes but the same exponents. The position and velocity distribution 
functions are only weakly anisotropic. 

PACS numbers: 05.70.Ln, 64.60.Cn, 64.60.My, 64.75.+g 
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I. INTRODUCTION 



There is on going interest in the growth kinetics of stripe forming systems. There has been 
progress via experimental and numerical studies of growth after a quench from an 

isotropic initial state. However the theoretical understanding of such systems remains limited. 
This is mostly due to the complexity of the defect structures generated during ordering in such 
systems. For example, in the Swift-Hohenberg model, there are grain boundaries, disclinations 
and dislocations generated in the ordering process. The co-existence of all these different defect 
structures has hindered the theoretical analysis of the striped phase ordering systems. In this paper, 
we study an anisotropic Swift-Hohenberg (SH) model, where only dislocations are produced in the 
ordering process. Our goal is to understand the statistical properties of these defects much as we 
now understand those properties for simple vortex producing models. 

There are formal arguments!^ that if we break the symmetry of the isotropic SH model by 
applying, for example an electric field, then the system can be mapped onto an an isotropic TDGL 
model. This suggests a L f^/^ growth law compared to much slower growth in the isotropic SH 
model. We find support for this hypothesis. 

Some previous studies have focused on the evolution of a few dislocations 
ll^ . Tesauro and Cross 0| studied the steady state climbing motion (move along the direction of 
stripes) of isolated dislocations both theoretically and numerically in several two-dimensional model 
systems including the SH model. They found that the wave number selected by dislocation climb is 
marginally stable only for potential models. Bodenschatz et al. 3| studied the climbing motion of 
dislocations with amplitude equations appropriate for systems with an axial anisotropy. The Peach- 
Kohler (PK) force ( the effective wave number mismatch) drives the dislocation motion, just as in 
Ref. M. Thev also consider the interaction between two dislocations together with the PK force. 
0„„:L al. & a a s..^^ t^e . „f . ne^a*. .a.eHal e.pe..enta,I. 



They introduced a gauge-field theoretical treatment to study the climbing of dislocations in a 
stressed background field where the PK force plays a role. The theory Q| predicts that climbing 
and gliding motions of a single dislocation are equivalent (after the proper scaling for the anisotropic 
system) and due to the PK mechanism. Braun and Steinberg studied the same experimental 
system. They measured the gliding motion of dislocations due to a pure interaction between the 
members of the pair without the PK mechanism. They found that the climb and gliding motion 
have different characters. 

Boyer [l^l simulated an anisotropic stripe forming model Q| based on the Swift-Hohenberg 
model. His model is more complicated than ours. In his model the stripes have two preferred 
directions and a zig-zag pattern is formed, and the dislocations tend to stay together to form large 
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domain walls. The author found that for small quenches the energy, the dislocation energy and the 
characteristic length normal to the stripes all scale as i^^/^ (+ for the characteristic length). He 
also found that for deep quenches the system was frozen. The pinning effect becomes important as 
the quench depth increases. The zig-zag pattern was experimentally realized in Ref. 

Here we study an ensemble of well separated dislocations in the context of domain growth. The 
motion of the dislocations in this model is highly anisotropic. They tend to move across the stripes. 
The average speeds across and along the stripes obey simple power laws in time with different 
amplitudes but approximately the same exponent. The distributions of the defect velocities along 
the two orthogonal directions have same form and large velocity power-law tails with approximately 
the same exponents. Two bulk measurements of the ordering, the decay of the effective energy and 
the number of dislocations, obey a simple power law in time with a logarithmic correction, as for 
the XY-model Q. 

The two dimensional isotropic Swift-Hohenberg (SH) model Q| is defined by a Langevin equa- 
tion 

^-^-«-). 

where ip is the ordering field, and the effective Hamiltonian is given by 



nm = I Sr \^-if + \ [(V^ + 1)^]' + i^^} , (2) 



where e is a positive constant. All the quantities in this paper have been put in dimensionless form. 
The noise ^ satisfies (^(x, t)^(x', t')) = 2T(5(x — x')(5(t — t'), where T is the temperature after the 
quench. In the following, we set T = which eliminates the noise term from the analysis. Starting 
from a random initial condition without long distance correlations, the SH equation Q generates 
stripes with period 27r. 

In the simulations for the isotropic SH model, we found P| that the grain boundaries' motion 
dominate the ordering dynamics of the system, which is different from what is seen in some ex- 
periments , where the disclination quadrapole annihilation is the dominant ordering process. 
Disclinations and dislocations are also present in the SH ordering system. In more recent experi- 
ments Ql on different diblock copolymer systems, defect configurations looking more like the SH 
simulations P| are found. The co-existence of different kinds of disordering defects makes it difficult 
to analyze simulations of the isotropic SH system. However, in an anisotropic SH system, where 
only point-like dislocations are present, the system should be easier to study. 

We make the SH system anisotropic by adding an additional term to the effective Hamiltonian 
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nm^nm^ \d'ri^\-^] , (3) 
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where 7 is a constant. The anisotropic term corresponds to applying an external magnetic field. In 
this case the SH equation now takes the form 

^ = .^ + (V^ + l)V-^3 + 70. (4) 

We studied the case where e = 0.1 and 7 = 0.4. Stripes generated by the above equation align 
along the y direction on the x-y plane. This configuration minimises the anisotropic term in Eq. 
iPJ. We start from a random initial condition for V'- After a very short transient time, the only 
defects left in the system are dislocations. Dislocation annihilation is the final ordering process. 

In Sec. II, we set up our numerical study. Then in Sec. Ill we study the time decay of the 
system energy. The stripe patterns and the motions of the dislocations are shown in Sec. IV. And 
we analyse the quantitative measurements in Sec. V. The speed distribution for the dislocations 
are shown in Sec. VI. 

II. NUMERICAL ALGORITHM 

We employed the usual Euler method to drive the system: 

i){t + M) = ^{t) + At (eV(t) - (1 + V^) + jdlil^it) - i^{tf) (5) 

We take in this case time step At = 0.02 and lattice spacing Ar = 7r/4. In the following sections, 
the numerical measurements are obtained from the system that is evolved by the Euler method. 

III. TIME DECAY OF THE AVERAGE ENERGY 

The first quantity we look at is a gross statistical measure of the ordering given by the average 
coarse-grained energy E = {7i)t- In Fig. 1 we plot AE = E — Eq, where Eq is the ordered value 
of E (known to be accurately given by Eq = — e^5/6, where S is the surface area of the system). 
These simulations were averaged over 528 runs. In agreement with the n = d = 2 TDGL model 
[itI . Q] we find a power law with a logarithmic correction characteristic of the annihilation of point 
defects. This is clearly consistent with a growth law exponent of z w 2. Although Fig. 1 does not 
give a good estimate for the exponent z, we conclude that z = 2 is the best value by taking into 
account Fig. 6 discussed below. Having established that the ordering is speeded up relative to the 
isotropic SH model, where z ^ 3, we can move on to look at the nature of the ordering patterns 
grown using this model. 
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Figure 1: AE v.s. time t after the quench. The energy AE is proportional to [i/ log(i/18.5)] 
averages are of 528 runs. 
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IV. DISLOCATIONS 



In Fig. 2 we show a typical ordering configurations of the anisotropic SH system at different 
times. Notice that the only defects produced are dislocations. 

We are interested in the path of each dislocation trajectory. These paths must be determined 



accurately enough such that we can compute dislocation speeds. In the case of the 0{2) TDGL 
model [ij, we were able to accurately determine the position of a vortex by finding the zeros 



(minima) in the order parameter amplitude. Here the situation is more complicated. As explained 
in reference P| the positions of defects in the SH model are located by maxima in the quantity 



A 



{Vc 



(6) 



where (p is the angle that the director h = VV'/|VV'| makes with some arbitrary direction and 
a is the index for different spatial directions. We showed numerically that if A > Aq then that 
site on the lattice can be associated with a defect. Here we need to determine the position of the 
dislocation with some accuracy. We have found that the expression 



(7) 



gives the position of the dislocation in the ath direction, and the sum is over all contiguous sites 
where Ai > Aq = 3.0. Using these procedures we obtain, for example, the set of dislocation 
trajectories shown in Fig. 3. The dislocations tend to move (glide) across the stripes and annihilate 
with each other. If we look more closely we see the oscillating behavior in the glide motion as shown 
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Figure 2: Typical configurations of an 512 x 512 anisotropic SH system at different times 
and top to bottom, the system was a.t t = 2500, 5000, 7500, 10000, 12500 and 15000. 
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Figure 3: The trajectories of the dislocations in a 512 x 512 anisotropic SH system. The dots are the 
positions where two dislocations annihilate. The stripes are along the y-direction as shown in Fig. 2. Most 
of the dislocations glide across the stripes. 



in Fig. 4 and 5. In Fig. 4 we also show all the sites that are associated with the dislocations. To 
obtain useful statistical data we will need to average over the different runs for our system. 

The oscillation in the glide motion of the dislocation is due to Peierls-like pinning forces. The 

n 

law of motion of a dislocation due to the Peierls-like force takes the form [U| 

_i dx 



dt 



f — p cos{kx) . 



Here v is the velocity across the stripes, /j, is the mobility. / is the external force per unit 
length applied onto the dislocation. In this context, / is caused by other dislocations especially the 
one that is going to annihilate with the dislocation of interest, p is the magnitude of the pinning 
force. The Peierls-like pinning force term oscillates with a period of 27r//c, which is exactly the 
stripe pattern period. We are interested in the interaction / between dislocations. In the next 
section, we will compute the the average v versus the separation distance between two annihilating 
dislocations. 

We plot the average number of dislocations Nd as a function of time in Fig. 6. Clearly it is fit 
w a power law with a log correction just as for AE and the isotropic TDGL result for n = d = 2 
13, Nd and AE' share approximately the same time dependence with growth law exponent 
z ^ 2. The energy per dislocation is almost a constant, which means the dislocations control the 
dynamics of the system. 
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Figure 4: The trajectories of two annihilating dislocations. The solid Hne denotes the center of the dislocation 
core determined by Eq. (7). The dots indicate the sites i which are part of the core used in Eq. (7). Thus 
this figure shows the motion of the dislocation core. 
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Figure 5: The position of a dislocation along the stripes versus the time after the quench. We use a 10th 
degree polynomial to fit the data and average out the oscillations along t/-direction. 
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Figure 6: The number of dislocations Nd for the anisotropic SH model as a function of time averaged over 
528 runs is given by the solid line. There is an excellent fit to the form {t / \og{t / a))~'' with a = 32.1 and 
b = 0.94. 

V. ANALYSIS 

To quantify the extent of cross stripe migration of dislocations, we measured the average distance 

in the climb and glide directions between two dislocations which are going to annihilate with 
each other at time to- In Fig- 7, we show v.s. t — to, where t is the time we measure the distance. 
We measured the components of average distance across and along the stripes. In Fig. 7, we can 
see that at any t — to the separation of two annihilating dislocations across the stripes is much larger 
than the distance along the stripes. This means the two dislocations tend to approach each other 
along the direction across the stripes. We notice that the average separation takes a power law form 
for the glide motion. The climbing motion as one approaches annihilation is more complicated. 

If we plot the average separation of defects heading toward annihilations versus time after quench 
we obtain the result shown in Fig. 8. Unlike in the NCOP TDGL case f deviates from a power-law 
behavior for long times and does not serve as a good measure of the growth law for the system. 

In Ref.|17| we found good scaling results and a reasonable theoretical model for describing the 
numerical results because there was a simple scaling relation between the average separation r of 
two annihilating defects and their relative speed u, u ~ r~^. The situation is more complicated 
here. In Fig. 9 we plot the average relative speed in a given direction versus separation in that 
direction. As the dislocations approach each other, their speeds increase. And when the distance 
between them is too small, our measurement is unable to follow the motions of the dislocations. 
So the data points in Fig. 9 with distances smaller than r = 16 are not reliable and should not be 
taken into account in the following analysis. 
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Figure 7: The average distance between two dislocations which are going to annihilate with each other versus 
time before the annihilation. The components across (glide) and along (climb) the stripes are measured 
respectively. 



Using only the reliable data in Fig. 9 we obtain 

Ua = Aar~'' (8) 

where b = 0.75 for both climb and glide, and Ag = 0.58, Ac = 0.06. The assumption Ua = Ua{ra) is 
only approximately correct. In fact Ua depends on both x-separation and y-separation distances, as 
is shown in Fig. 10a. Correspondingly, there is a strong correlation as shown in Fig. 10b between 
the glide separation dg and the separation distance d. The climb distance is correlated with d for 
small enough separations. 

If we plot the average relative speed versus time to annihilation, Fig. 11, we obtain the approx- 
imate power-law results 

Ua{t)-{to-t)-'/' , (9) 
which is consistent with z = 2 m both directions. 
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Figure 8: The average distance between two dislocations which are going to annihilate versus the time after 
the quench. The components across (glide) and along (climb) the stripes are measured respectively. 

VI. PROBABILITY DISTRIBUTIONS 

We can next turn to the associated probability distributions. The first is the separation prob- 
ability distribution. This is the probability P^it) that at time t two annihilating dislocations are 
separated by a distance r^, where a corresponds to climb or glide. We assume that P^{t) takes a 
scaling form if we plot it versus r a/ fait) as shown in Fig. 12, where the average fa is shown in 
Fig. 8. We obtain a reasonable scaling form but unlike in Ref. 0| we do not find an algebraic 
large separation tail. Instead the scaling function appears to decay exponentially. Notice that the 
scaling forms are roughly independent of direction. We fit the separation distribution function with 
y = aQX^/{l + 012;^+^)'^, where ao, oi, b and c are parameters. This function is the most general 
form for the separation distribution obtained from Since the fitting function has power-law 

decay tail, which is different from the exponential decay of the real data, we are unable to fit the 
large separation tail. When we fit the data, we restrict the value of c to be between 1 and 5. The 
resulting value of c is 4.9. If we extend the upper limit on c, the best fit for c increases. However 
we found that b is always close to 0.7 whatever the upper limit of c. 

We turn next to the statistics governing the relative speeds. The average speeds are shown in 
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Figure 9: The average relative speed between two dislocations which are going to annihilate versus the 
separation between them in that direction. The best fit for each set of the glide data is given by u(ra) = 
Aar~^ with ha = 0.75 approximately for both glide and climb. 

Fig. 13a to be given approximately by t^^/^ for both directions in agreement with the scaling ideas. 
If we fit the climb trajectories in a way that averages over the oscillations in the climb component 
we obtain the average climb speed which is in better agreement with the result, as is shown 

in Fig. 13b. 

Finally we plot the speed distribution function in Fig. 14a. Clearly we have large speed power- 
law tail. We find roughly that both glide and climb motions have a f ~^ large speed power-law tail 
as shown in Fig. 14b after we average out the oscillations on the y-component (along the stripes). 
The distribution functions are sensitive to how we treat the oscillations in the climb data. This 
may account for the higher tail exponent 3.9 shown in Fig. 14a for the climb data. We also use a 
quite general form for the speed distribution function obtained from Ref.jl3| to fit the data. The 
function is y{x) = aQX~'^~^^'^^'^~^'^'^~^^/^ / {1 + aix^'^'^'^'^^^Y . We require that y(0) to be nonzero. So we 
must have —2 -|- (c(6 -|- 1) — l)/6 = 0. The parameter b is the same b in Fig. 9, with a value 0.8. 
So we have c = {1 + 26)/(l -|- 6) = 1.44. So we fix the values of c = 1.44 and b = 0.8 in the curve 
fitting. 

Again the scaling forms in the two directions are, to within our accuracy, equivalent. 
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Figure 10: (a) The speed Ux and Uy depend on both x and y. (b) d is the distance between two annihilating 
dislocations, dg is the separation distance across the stripes (glide motion), dc is the separation distance 
along the stripes (climb motion). And d^ — d^ + dg. 

VII. CONCLUSION 



The kinetics of the anisotropic SH model are conceptually simpler than for the isotropic SH case 
since there is only one disordering defect. In the simplest picture one has a set of point defects 
ordering in a fashion similar to a collection of vortices in an XY model Q| but with anisotropic 
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Figure 11: The average relative speed between two dislocations which are going to annihilate versus the 
time before the annihilation in the glide and climb directions. 




Figure 12: The separation probability distribution. The tail of the distribution has an exponential form, as 
is shown in the insert. We fit the separation distribution function with y = a^x^ /{I + aix^'^^Y , where oq, 
ai, h and c are parameters whose values are given in the figure. 
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Figure 13: (a) The averages of the speeds, the transverse component (across the stripes) and the longitudinal 
component (along the stripes). They all obey simple power laws Va ~ t^" ■ In (b), we use 10th degree 
polynomial to fit the trajectories of cHmb motion and average out the oscillations on climb motion. The 
obtained power-law exponent is closer to 0.5 in this case. 
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Figure 14: The distribution of the defects' speed for glide and climb directions. After scaling the data at 
different times (from 300 to 8000), we collapse all the speed data to one curve. The motions of gliding 
and climbing are measured separately. In (b) , the oscillations of the climbing motion along the stripe are 
averaged out. We also use the 10th degree polynomial fit to smooth the ghde motion. 

scaling. One has a scahng length La{t) ~ A^t^ and power-law behaviors for the decreasing energy, 
number of defects, and average defect velocity with x ~ 0.5. 

The picture of annihilating point defects with a growth law of t^/^ is roughly true for our 
anisotropic system and in this sense the system is similar to the 2D XY model with the annihilation 
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of vortices. Both have bulk properties AE and N^i which have the same ordering time dependence. 
Both systems have large velocity power-law tails which show behavior. One simple result is 
that the scaling functions P{r) and P{v) are nearly isotropic as shown in Fig. 12 and 14. This 
is true despite the fact that the motion is highly anisotropic. It appears that the annihilating 
defects organize themselves in such a way that their relative velocity is radial. Once this is true 
it is not important what the angle is between the relative velocity and the stripes. As one looks 
closer the analogy breaks down. The separation probability distribution shows a clear power-law 
tail in the XY model but not in the anisotropic model. The average distance between annihilation 
defects in the XY model serves as a good measure of the growth law f ~ t^/^. This is not true in 
the anisotropic model where f ~ t^-^"^ . One interpretation of our results is that the hypothesis of 
independent pairs of dislocations breaks down at a much greater distance than for the XY model. It 
is also possible that the independent pair mechanism used in Ref . \r\ works less well in this system. 
There may be correlations among different pairs of dislocations. In any event our results can not 
be simply explained by a rescaling of the glide and climb directions. 
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